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Abstract 

A irrotational solution is derived for the steady-state Navier-Stokes 
equations that approximately satisfies the boundary conditions for flow 
over a finite flat plate. The nature of the flow differs substantially from 
boundary layer flow, with severe numerical difficulties in some regions. 



An analytic function having the form 

f(z) =v + iu, (1) 

leads to an exact solution of the two-dimensional steady-state Navier-Stokes 
equations, i.e., 

p u ■ Vm = — Vp + /iV 2 «; 

V-it = 0. (2) 

This occurs because (1) leads to a velocity field having the properties 

V 2 « = 0; (3) 

and 

V-m = 0; 
V x u = 0. (4) 

This is a subset of the generalized Beltrami flows 1 , and note that (4) are 
the Cauchy-Riemann equations for (1). Substituting (3) and (4) into (2), and 
using the identity 

u-V u — (V x it) x it + ^V-u 2 (5) 
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leads to 
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2 Vu 2 = -Vp, (6) 



or 

,2 



The main difficulty with (1) involves finding a function f(z) satisfying useful 
no-slip boundary conditions. Consider the function 



sin (2ir (z — zq) / L) idxo 
A+[B sin (2tt (z - z ) /L) - A] e (^ z o) 2 A 



(8) 



Here z = x + it/, z = x + *yo and y = 0. When y — > oo 

sin (27T (z — zq) / i) idxo 



Making use of the identity 



= 0. (9) 



sin (2tt (z — z ) / L) = cos (2ttx /L) sin (2nz/L) (10) 

— sin (27rx / L) cos (27rz /i) , 

the terms cos(27rz/L) and sin(27rz/L) become large as y — > oo; however, 
sin (27ra;o / i) and cos (2nxo / L) are precisely zero when integrated over < x < 
L. 

When y = 0, 

= lim I / L sin(2 7 r(,-.x )/£)^ 

' e-^oLJo A+ [B sin (2tt (.t - a; ) /L) - A] e (^-^o) 2 /^ 

When x ^ xq, the denominator in (11) diverges. However, as x — > xo, the 
denominator approaches Bsin(27r(x — x ) /L), so that as e — > 0, u approaches 
a rectangle function between x = and x = L and t> approaches zero, meeting 
the boundary conditions for flow over a flat plate. 

Although the integral (8) can be difficult to evaluate in some regions, soft- 
ware packages seem to have less trouble converging when B >> A, at least for 
y = 0. Choosing e = 0.0001, A ^ 5.6 x 10~ 4 , B = 1Q0A and L = 1 m satisfies 
the flat plate boundary condition as shown in figure 1. As e — > 0, u becomes 
more step- like at x = and x = 1. Numerical approaches require e to be finite, 
leading to regions at the plate edges where the boundary conditions are not 
satisfied. These regions can be made as small as required by reducing e. 

The boundary layer assumptions, i.e., >> 1^ and = are useful for 
analyzing flow over a flat plate, but they break down in the region at the leading 
edge 2 , when U^x/v < 10000. This region of non-validity can be used to define 
an acceptable region where the flat plate boundary conditions are only satisfied 
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Figure 1: The u approximation to the flat plate boundary condition at y = 
for e = 0.0001. 



approximately by (11) and thus determines a practical value for e. When 
U = lm/s and v = 10~ 6 m 2 / s, the extent of non-validity is approximately 
defined by x < 0.01m. Setting e = 0.00001, A = 5.9 x 10~ 5 , and B = 100,4 
leads to the boundary conditions being satisfied to within a tolerance on the 
order of 10~ 6 at the edges of the regions -0.01 < x < 0.01 and 0.99 < x < 1.01 
as shown in figure 2. 

Although the flow described by (8) exactly satisfies the boundary conditions 
for a finite flat plate when e — > (and approximately for small but finite e), the 
flow field differs substantially from boundary layer flow. Figure 3 shows the re- 
gions having nonzero velocity components. As y increases, the integrand evolves 
into a sinusoidal-like function, until it becomes sufficiently sinusoidal that u and 
v both approach zero. Just below the upper "zero" line, small perturbations 
in the sinusoidal-like integrand functions lead to nonzero velocities. In some 
regions, either the first term or the second in the denominator in (8) dominate, 
allowing the other to be neglected. The latter is true when y << 1 above the 
plate, allowing the velocity field to be accurately approximated. When nei- 
ther terms can be neglected, evaluation of the integral becomes tedious because 
of very large integrand magnitudes, leading software packages to either fail to 
converge or converge to incorrect solutions. Numerical determination of the 
velocity field for many regions remains a challenge. 
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Figure 2: The u approximation to the flat plate boundary condition at y = 
for e = 0.00001. 
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Figure 3: Upper and lower "zero" lines. All nonzero velocity components are 
confined to the region between these lines. 
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